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1. Introduction 

The Standard Model(SM) of elementary particle physics has been proved to be a remarkable 
successful description of high energy physics phenomena.However, the origin of neutrino 
mass is still an unknown fundamental problem[|3l 01. In SM neutrinos are being introduced 
as a zero-mass.The discovery of flavor conversion of solar and atmospheric neutrinos[|5l |6l 
has established that neutrinos have nonzero mass and they mix among themselves, therefore 
providing the first evidence of new physics beyond the standard model. 

On the other hand, there has been growing interest in the anomalous magnetic moments 
of during the past few years. The most recent determination of a^ in the SM is[|71 



The difference between them is 3(7. Though 3<7 can not be regarded as the powerful proof of 
the new physics,it is possible that the difference will become larger with the development of 
experiments. So, we maybe get the proof of the new physics later. 

A model with the above desired feature built in has been suggested recently by Hung 
|l9l. The model keeps the SM gauge group albeit in a 'vector- like' manner: Mirror fermions 
charged under SM gauge group, of which right-handed neutrinos are a member, and two Higgs 
triplets in a way that preserves the p parameter to unity, plus a Higgs singlet. In particular, a 
right-handed neutrino which is sterile and are required to have a mass now becomes a member 
of a weak doublet of mirror leptons. A tiny Dirac mass for neutrinos is offered by a scalar 
singlet whose vacuum expectation value is not necessarily associated with the electroweak 
scale, while a Majorana mass of order the electroweak scale is provided by a scalar triplet. It 
is conceivable that these new leptons could be discovered at high energy colliders in the near 
future and the rich lepton flavor structure could be observed in low energy processes. The 
weak charged couplings are generally non-unitary with or without restricting to the subspace 
of light leptons, and flavor changing neutral currents (FCNC) occur in a way that is controlled 
by the weak charged couplings. 

The paper is organized as follows. First, we will introduce the Hung model briefly; for 
a full account of the model, we refer to refs Then we show in section 3 how to calculate 
the anomalous magnetic moment of jU in the virtual transition /i — )• fiy in gauge, and 
demonstrate the final result. Our result is summarized in the last section. 

2. Hung's Model 

We start with a brief description of the model relevant to our later analysis; for a full account 
of it, see Ref.[91. Considering three generations,the SM and mirror leptons with quantum 
numbers under the gauge group SU (2) xU{1)y are: 



af^ = (11659183.4±4.9) x 10"^" 
The total SM prediction of ji differs from the experimental valueflU 
af'P = (11659208.0 ±6) x 10"^^. 



(1) 



(2) 
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where the subscripts L, R refer to chirality and the superscript M to mirror, and the first 
number in parentheses stands for the dimension of representation under the gauge group 
SU (2) . Because of anomaly cancelation,the quark sector also has mirror partners which are 
of no interest here.Besides the SM scalar doublet 4>, the model contains the new scalars 

0(1,0), z (3,2), (4) 

plus an additional triplet ^ (3,0) that together with x preserves the custodial symmetry [[TOl 
but is irrelevant here. 

The Yukawa couplings of leptons are, with the generation indices suppressed, 

- = +yMW^fL + h-c, 

-^^ =x^JZF/0+x/^/f0+h.c., 

(5) 

where Xj/^ = ^/'i//*, = i"/^Y^, and 



-^x =lzMiFi'niT^)xFi' +h.c.. 



1 



V2 [ V2x^ -X 
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Note that a potential Majorana coupling of x to Ft is forbidden by imposing an appropriate 
f/(l) symmetry . Suppose the VEVs have the structure: 




(<^>)=vi, (Z) 




(7) 



where V2,3 contribute to the masses of weak gauge bosons and are naturally of order the 
electroweak scale while vi is not necessarily related to it. In the basis of /, /^,the charged 
lepton mass terms are 



m 



f 



V 

while the neutrino mass terms are 

1 
2 




( ^ 



+ h.c.. 



(8) 
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c/>n 
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+ h.c.. 



m„ 



VlXf 



VlXp V3ZM 

We denote the mass eigenstate fields of the charged and neutral leptons by ij. 



(9) 



V; 



respectively, with j = 1,2,3 for those that are mostly ordinary leptons and j = 4,5,6 for 
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those that are mostly mirror. The seesaw mechanism operates for a Majorana mass of order 
the electroweak scale and Dirac proportional to Vi which can be chosen small.This relaxes in 
some sense the tension in ordinary seesaw models between the generation of light neutrino 
mass and the observability of heavy neutrinos at colliders. The mass matrices are diagonalized 
by unitary transformations as follows 




(10) 



where the unitary matrices Xl^r and Y satisfy 

xlmfXR = m£ = diag(m^.) , F^m„F = my = diag(mv,) , (11) 

with the mass eigenvalues m^. y, being real and nonnegative. There is a constraint on their 
masses from the zero texture, 

6 

Y,my,YikYjk = 0, for/,; = 1,2,3. (12) 

k=i 

Splitting each unitary matrix into two blocks that contain the upper and lower rows 
respectively, 

the following relations will be required later: 

fL,R = XIrIl^r, ff^R = V, HL = 7"* Vi, nf = yW. (14) 
where the unitarity constraints become 

^a^a* = XpC^ = I3, XpC^^ = O3 (for a = L,R), 

yuyu+ ^ ydyd+ ^ yUj.d+ ^ Qg, (15) 

and the zero texture constraint only acts on Y^: 

FVyy"'^ = 03. (16) 

The above diagonalizing matrices will enter the gauge interactions of leptons.Some 
algebra yields, 

= 82 {jw^W+ + j^^W- +J^Z^) +eJ^mA^, (17) 
where the currents are (Pl,r = 1 =F 75/2) 

V2j+^ = VY^ {VlPl + VrPr)£, 

cw4 = {VlVIPl + VrVIPr^ V - ^-If {vIVlPl + vIVrP^ i 

JL^-ly^i, (18) 
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with V = v^; + Vl = v^; + = Vl + V^,and cw = cos Bw,sw = sinSvK with 9^ being the 
Weinberg angle. To relate the matrices Vl^V/? to Xa {a = L,R),Y ,\i is convenient to decompose 
the latter into the up and down 3x6 blocks, 

then 

VL = Y''''XlyR = Y^'^Xl (20) 

with VIVr = 0. These matrices are generally non-unitary and the deviation from unitarity 
induces FCNC in both sectors of neutrinos and charged leptons: 

Viyl = Y^TyU*^ y^yl ^ ydfyd^ 

vIVl = Xfxl, V^Vr = (21) 

The matrix blocks X^, X^ do not enter charged current interactions since are SU (2) 

singlets. It is important to note that the matrices in charged currents, Vl, Vr, are generally not 
unitary and there are flavor changing neutral currents arising from the nonunitarity. 

If we go on with our calculations, we will also need the Yukawa couplings of the would- 
be Goldstone bosons (GB's). Although the original scalar fields mix in a complicated manner 
via the terms in the scalar potential, GB's are independent of mixing details due to gauge 
symmetry and their Yukawa couplings only involve matrices that already appear in gauge 
interactions. We list the Feynman rules as follows: 

VAG+ : [miVtaPL + VUPR) -maiytaPR + Vi'M] , 

y2mw 

laViG- : [m^t^PR + Vl'^'PL) -ma{Vt^ Pl + V^'Pr)] • (22) 

We give below other Feynman rules that will be required later (all momenta incoming): 

i^e{-i){k+-k-)^ = ie{k+-k-)^, 
iemwg^v, 

TajS/x = gaj5{k--k+)^+gp^{k-k-)a 
+ g^a{k+-k)p. (23) 



A^G+ik+)G-{k^) •-•2- 



3. Calculation in gauge 

Now we begin to calculate the anomalous magnetic moment for ji in both W diagram and 
Z diagram. As discussed in [|2l, there are two types of contributions at the one loop level that 
are mediated respectively by charged current interactions of and flavor changing neutral 
current interactions of Z^, corresponding to unitarity gauge and gauge.For the charged 
current diagrams we find that they involve a triple gauge coupling which is more divergent 
in the ultraviolet and the flavor mixing matrices in the charged current are not unitary.Thus 
in such a circumstance, it is highly desired that the calculation should be done in a safer R^ 
gauge whose ^ dependence is canceled as expected. 
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Figure 1. Diagrams that contribute to the decay jx ^ jxy. 



3.1. The contribution of W diagram 



First of all,we should calculate the amplitude for the decay — ll{p — q)y{q,£) in W 
diagram whose Feynman diagrams are shown in Fig l.The anomalous magnetic moment of jU 
is: 

z _ ^-2 



^2(^^ = 0). 



where 



2m,, 



(24) 



(25) 



The Lorentz and gauge symmetries dictate that the amplitude has the following structure [[TT|: 

r^u^iOiyq''[x + yY5]£^u^, (26) 

so that we can concentrate on the iOx^q^ terms to pick up the coefficients, x and y. 

We shall keep only terms up to the linear order in the muon mass, m^. These are indeed 
very good approximations. The four diagrams give the following on- shell amplitudes ;for a 
full presentation of our calculation details ,see Appendix: 

/ 1 



T{a) = eg2Uf^{io^vq^] 



+miB 



{An 



>2 ,^2 



m 



w 



m^^A 



M{a)+P{a)ln[n] + Q{a)\n[Q 



M'{a)+P'{a)\n[n] + Q'{a)\n[^ 
+m^B\M"{a)+P"{a)\n[ri 
+m^D iMoia) +PD{a) ln[r/] + Qoia) ln[<^ 



(27) 



T{b) = eglu^iio^vq 



{An 



i 1 

?2'Z2 



W 



m^A[M{b)+P{b) Inn + Q{b) In^] 

+miB [M\b)+P'{b) Inn + Q\b) In^] 
+m^D [Md (b) +PD{b)\nn + QDib)ln^]}u^, 



(28) 



T{c) = eglu^iio^vq"^ 



i 1 



{4n 



,2 ^2 



m 



w 



m^A 



M{c)+Pic)\nn + Qic)\n^ 



+miB 



M'{c)+P\c)\nn + Q\c)\n^ 



+ m^D 



MD{c)+PD{c)\nn 
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Figure 2. Diagrams that contribute to the decay — > /xy. 
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(29) 



T{d) ^ eglu^iio^vq"^] 



(An) 



i 1 , 



m 



w 



M{d)+P{d)\nri + Q{d)\n^ 



+miB 
+m^D 



M'{d) +P'{d) Inn + Q\d) In ^ 

Mo ( J) + Pd (d) In n + 2/) ( J) In ^ ] } . (30) 

The loop functions such as M{a),M{b),M{c),M{d).... and other functions are listed in 
Appendix. 

Collecting the contributions from the four graphs yields the final answer: 
Aw = I,iVl,^Vun+V;,^VR,^) = Dw, 



where ^ = .The loop functions are 

1 



V2 8ml 



6{i-n) 
1 

1 
1 



10 - 43r,- + ISrf - 49rf + 4rf + 18rf In n) 



{-4+15r,-12rf + + 6r} In r,) , 
[I3r2 - lAn + 11+ 2r,(5 - 6r,) Inr,-] , 
[7 - 34r,- + 33r2 - 1 Orf + 4rf -l%r} In r,-] . 



6(l-r,)4 

It is good to see that the t, dependence is completely cancelled as expected. 



(31) 



(32) 
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3.2. The contribution ofZ diagram 

We begin to calculate the amplitude of diagram Z whose Feymnan diagrams are shown in Fig 
2.We can use the same method in the diagram W.Thus we can give the contributions from the 
two graphs in Z diagram;for a full presentation of our calculation details ,see Appendix: 

2 

rmz] = 4 J^^2 (4^)2 vg^ ) {^nA' [MM{b) + PP{b) In r« + QQ{b) In ^ ] 

+maB' [MM'{b) +PP'{b) Inra + QQ'{b) ln<^] 

+m^D'[MMo{b)+PPo{b)\nra + QQD{b)ln^]]u^, (33) 

2 

T[bGo] = ^^j^u^,{ia^,yq'')[m^,A'[MM{c)+PP{c)\nra + QQic)\n^] 

+maB' [MM'{c) +PP'{c) Inra + QQ'{c) ln<§] 

+m^Z)' [MMd (c) + PPd (c) In ra + eei) (c) In ^ ] } Mju , (34) 

Using (31) and (32), we finally get the contribution of the diagram Z to the anomalous 
magnetic moment of ju: 

2 y/lGfCw 



a 



mlA'z^z{ra)+m^maB'z^z{ra) +mlB'zJ^z{ra)+mlD'z^z{ra)\ , 

a 
a 

where 

4(1 -ra) 

^^z = — [-4-ra + 8r2-3r3 -6r2lnra], 

2(1 -ra)'' 

^z = 77-^ - 13ra + 1 1 + 6ra(5 - 6ra) Inr^] , 
^z = ^^^^^^^^4 [- 14 - 6r« + 57r2 - 54r3 + 1 Br^ - In r«] . (36) 
We find the result is also nothing to do with ^ as expected. 



4. Numerical analysis 

We have calculated the W and Z diagrams in gauge.An interesting technical point is in 
order. We can also work in unitarity gauge and it should be simplest.But in the W diagram 
there is more ultraviolet divergent due to the triple gauge coupling. Thus there is no guarantee 
in this case that the order of removing the ultraviolet regulator commutes with that of taking 
the unitarity gauge limit.As a matter of fact, although the diagram is convergent in both 
unitarity and gauges, there is a finite difference in the terms linear in the lepton masses 
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Table 1. The anomalous magnetic moment of jU(a|j x 10~") is increasing with the increase 
of the heavy neutrino mass when the unknown a mass is a fixed value at first, while it is 
decreasing when the mass of the heavy neutrino is over 700 GeV. 
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4.58481 
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300 


5.68684 
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5.69074 


5.69079 


5.69080 


5.69079 


5.69078 


400 


6.94412 


6.94722 


6.94784 


6.94802 


6.94807 


6.94807 


6.94807 


6.94805 


500 


8.26680 


8.26990 


8.27052 


8.27070 


8.27075 


8.27075 


8.27075 


8.27073 


600 


9.62381 


9.62691 


9.62753 


9.62771 


9.62776 


9.62777 


9.62776 


9.62775 


700 


11.0017 


11.0048 


11.0054 


11.0056 


11.0056 


11.0056 


11.0056 


11.0056 


800 


12.3934 


12.3965 


12.3971 


12.3973 


12.3973 


12.3974 


12.3973 


12.3973 


900 


13.7950 


13.7981 


13.7987 


13.7989 


13.7989 


13.7989 


13.7989 


13.7989 



between the results obtained in the two gauges. This caveat is restricted to the mentioned 
terms because terms of a higher order are convergent enough to allow the free interchange 
of taking the Hmits. Considering this,we should work in gauge make sure not to have any 
problems. 

Using (29) and (3 3), we finally obtain the anomalous magnetic moment of jU: 

a^ = aj^al = ^ = F^iq^ = 0). (37) 

The above anomalous magnetic moment involves several neutrino masses,the unknown a 
mass and many mixing matrix elements.In our later numerical analysis, we shall make some 
approximations. First,the light neutrinos can be safely treated as massless in the diagram 
W. Then, ^ ^ = ^ n,^w = ^ term ^ can be dropped 

because of m,-. Second, for the heavy neutrino mass we choose m/, = 200,300 up to 900 
GeV.We use mw = 80.2 GeV.Thus n = 6.21887.Then,^n/ = ^ 0.0883211,^,^ = ^ 
-0.350906„^VK ^^w^ -0.320226,^vi/ ^ ^ 0.920466. 

As a bonus of the approximations, the anomalous magnetic moment of ju in the diagram 
W depends on the products of matrix elements summed over light and heavy neutrinos 
respectively: 



W = Y.^yl)e0L)iv.. vi = Y.^Vl)e0R)iv.. (38) 

and similarly for V^^ 3 4 with / summed over 4,5,6. 

We find an algebraically simple case after some inspection. Suppose the upper-right 
3x3 block of Y is real. In this scenario , our special neutrino spectrum (three almost massless 
plus three almost degenerate and heavy) implies that the two off-diagonal 3x3 blocks of 
Y vanish, the lower-right block is trivially identity and the upper-left one is unitary. Then, 
W ~ (•*l-*z.)eM' ^1 ~ ~^^j^R)e\jL while all others vanish, where xl^ are the upper-left 3x3 
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blocks of Xi R respectively. Since we have no idea of their magnitudes, we sample randomly 
the real and imaginary parts of V/, Vj* between —2 x 10^^ and +2 x 10^^. 

Thus for numerical analysis we can assume Aw,Dw ~ 10^^, ~ 10^^ with i summed 
over 1,2,3 while By^/ — with i summed over 4,5,6. We also use a — l/ 131 .OA, m^, — 0. 1056 
GeV. 

We proceed to the diagram Z.First ,we can also assume the unknown a mass ma — 
200,300 up to 900 GeV.Then,ra = 4.81048.We use mz = 91.187621 GeV from PDG.Thus 
,we get :^z ^ -3.15976,^z ^ 2.28058,^z ^ 6.91432,^z ^ 1.90094.We also make some 
approximations: 

(4;^aM)<< ^ 2(v2v^)^^(y2y^)«;, -4sl{VlVL)a^ ~ 2(1 -24)10-^ 

Thus we obtain the different value of the anomalous magnetic moment of jU which has been 
shown in Table 1. 

It is interesting that with the fixed unknown ma the anomalous magnetic moment of jU 
is increasing with the increase of the heavy neutrino mass at first, while it is decreasing when 
the mass of the heavy neutrino is over 700 GeV. 

5. Conclusion 

The subject of the anomalous magnetic moments for the muon is an exciting and fascinating 
topic because it represents the best compromise between sensitivity to new degrees of 
freedom describing physics beyond the standard model and experimental feasibility. We 
have calculated in gauge the anomalous magnetic moment of jU in a model suggested 
recently. Because of the rich flavor structure of the model, the weak charged currents involve 
nonunitary mixings between the neutral and changed leptons and contain both left-handed and 
right-handed chiralities.lt is conceivable that these new leptons could be discovered at high 
energy colliders in the near future, while the rich lepton flavor structure could be observed in 
low energy processes. 
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Appendix A. The calculation of the contribution from W diagram 

We start with diagram (b): 



11 



where 



[k^-^ml][{k + qY-ml 



(A.1) 



N^, = Ya iyl)^iPL+{yl)^,iPR 



1 



Iji+^-mi 



:Um-(i-0 



{k + q)^ik + qy 
{k-\-qY — ^m- 



(A.2) 



To simplify A^^, we apply u^4 = u^{—mf^ i-c, we can make the replacement, ^ —)> ^ — m^, 
when 4 is on the far left. Note that the term q^x has no contribution.Now can be split up as 
follows : 



iV„ =1 + 11 + 111 + VI, 



(A.3) 



where 



I = 



n 



ni = 



{Vl)^iPR+{Vl)t,iPL 

-{\-^)[{yl)t,iPR + {vl)^iPL 
-(l-0(V;P/e + V;PL)m, 



^2 _j_ p2 _ yyj2 



2' 



ik + qY-Ww' 
k^{J^+^ + mi 



lW={\-^)iylPL + vlPR)m^- 



k^{l^+^ + mi) 



[{k + pY-m}] [{k + qY-^ml] 
These three terms will be calculated separately below. 

The term I gives the following contribution to the loop integral: 

{vl)^iPR + {vl)^iPL] Y^ill+^ + nii) 



(A.4) 



d'^k 



X 



(A.5) 



{2kY [k^ + p2 _ ^2] [^2 _ ^^2^] [(^ + ^)2 _ ^2 

mi [iVL)i^PL + {VR)i^PR^ - m^ (^{VL)i^,PR + {VR)i^PL) 

As mentioned before, we only need pick up the {iOp^yq^) terms. Using ^u^x = m^u^i, we can 
make the replacement / — )■ when ^ is on the far right. Since we work up to the linear 
order in m^, it is sufficient to expand to 0{p^). The expansion yields several types of terms. 
By symmetric integration, we get y^l^^—lk- p) — )■ —\y^]/)k'^ which has no contribution to the 
desired Lorentz structure. Similarly, the other two terms are, y^l^{—lk-q) —^7ix4^'^ 
j{io^v<l^)k^- The contribution of the I term is summarized as follows: 

k^ 



{b)i = ^eglu^,{iG^,yq')miB I — ' 

[k^ - mf] [k^ - t,m^] [k^ - m^] 



-^eS2Uti{i(y^vq^)m^A 



k[k^-mj][k2-^m^^][k^-m^^]^ 



(A.6) 
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where 

f f d^k 
Jk^ J (2^' 

i 

5 = E {{vl)^i{VR)i^PR + {vi)^,i{VL)i^PL) . (A.7) 

i 

We find that the term II has no contribution to (b). 
In calculating the term HI, we will use: 

kull^ik- p){k-q) 
1 1 

k^^{k-p){k-p) 

^ Y^{ia^^q'')]/>k\ (A.8) 

Note that adding the subscript * to means Pr, Pr — ^ Pi when Pj^ and Pr are on the far 
left. Thus we get the contribution of the term 111 to the loop integral: 

(^)m = ^eglu^{ia^vq^)mfm^A{l 



{k^-m}Y{k^-^ml)\k^-ml,) 
1 



(^2 -m2)2(fc2_ ^^2^)2(^2 _ ^2^)2 
1 



{k^-m]f{k^-^ml)\k^-m^^) 

+ \e8lu, M)mMl - ^ ) / (,2_^2)2(,24l2,)2(,2_^2^) 

+eg2U^i{io^vq^){\ - ^)m]m^D / 

Jk 



(V III 

J_ k^ 

24 (^2 _^2)2(^2_ ^^2^)3(^2 _^2^) 



\_ k^ 1 

^24 {k^ mj)^k2 - Ww)Hk2 - mi)^\ ^^'^^ 
We can use the same method in the term IV and we get the contribution of it: 

(^)iv = esWiia,,q^)il -^)mjm,Di-l) -^m\)\k^ -m^i "'- ^^'^'^ 

where 

^ = E (^i)MK^i?)iM^/f) • (A.11) 

i 

Now we can get the contribution of diagram b to the loop integral: 

T{h) = eglu^{io^vq^){m^AF+miBG + m^DW)u^, (A.12) 
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1 a 11 



4^-1 24<^-l 

, 1 n 



r-l 
1 1 



h 



1 n 



1 



24<^-l V<^-1 

2 1/4 + ^ 



2/3 



' 24^-1 n-1 7" 12(r,--l)(^-l) ' 12 

r ^ ^ 1 ^ ^ 1 1 ^ 1 ^ 2., 



^7 + — r,m^(^/9 



1 



24(^-l)(r,-l) ^-1 



1 



^-1 



1 

"24^24(<^-l) 
1 



24 24(^-1) 24(r,--l), 



^4 + ^m^n^Ig - ^m^n^^ho. 



(A.13) 



Note that = /m^Jn are integral functions listed in Appendix C. Substituting the /„ 
functions into T{b) gives: 



T{b) 



{An) 



'—-Legla^(i(y^^qy)Im^A [M{b) +P(&)lnr; + (2(&)ln<^] 



m 



(A.14) 



+miB [M'{b) + Inr,- + Q!{b) In^] 
+mjuD [Mo (&) + P£, (&) In r; + Qi, (&) In ] } , 

where M(Z7) , /'(Z?) , ^(Z?) , . . . are listed in Appendix E. 

Now we proceed to diagram (c) and (J). Following the steps which are entirely similar to 
those in the calculation of diagram [b), we can obtain the contributions to the amplitude from 
them: 

/ 1 



T{c)=eg2U^{i0^yq 



+miB 



{47t) 



2 ^2 



m 



w 



M{c)+P{c)lnri + Q(c)ln^\, 



M'{c) + P'{c)lnri + Q'{c)ln^\ +m^D Md{c) +PD{c)lnr, 



+Go(c)ln<^ 



(A. 15) 



T (d) = eglu^ {io^vq") \ \m^A \M{d) + P{d) Inr,- + Q{d) In ^ 

{An) ^ '- 

+miB \M'{d) + P'{d) Inr, + Q!{d) In t, 
+m^iD \Md (d) + PDid) In n + Qd (d) In ^ 1 1 , 



(A.16) 

where M(c) , P(c) ,M{d) ,P{d), ... are Usted in Appendix F and G. 

We finally come to the diagram (a) which is the most complicated one due to the 
appearance of a triple gauge coupling and double gauge boson propagators: 

- -^{a) = u^j ^^Ya [{Vl)^,iPL + {Vi)^iPR\ uA_^ yp [{yL)it,PL + {VR)i^PR] ^ap^ 
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X 



— m 



m 



w 



{k + qy-m^ 



X 



{k + qY-^w? 



W 



(A. 17) 



The above can be split into three terms corresponding to the product of the two propagators: 



— {a) = {gg) + {k) + {k + q) 
^82 



(A.18) 



where 



(gg) = 



d'^k 

X 



{vl)^,iPR+{vl)^,iPL\/ 



]jl+^-m, 



Y'[{VL)inPL+iVR)i^PR] 



a/3/1 



[k^-ml][{k + q)^-m 



(A.19) 



1-^ 



k'T, 



X 



aPu 



[A:2 - ml,] [k^ - ^mlr] [{k + qf - m 



(A.20) 



The fourth term from the product of propagators has been discarded since it does not 
contribute to the on-shell amplitude. 

We start from the apparently easiest (actually the most complicated) term (gg). Simplify 
by aiming at the leading terms linear in m^u : 

Note that the term q^j^ and the term can be dropped because they have no contribution to 
{gg) and we use 7(^(^+^)7" = — 2(^+^).Thus we get: 



Ijt+p-mi 



{k + pY — rm 



Alik^ + Ik^ {24 - 3mi) + 7^^(2^ -P) + y^i24{-^ + mi) + m,-(-^)7^ 



(A.23) 



Substitute (A.23) into (A.19) and again we should pick out (iOp^vq^) terms. By using the 
same method as in calculating {b)i we get the following contribution to the loop integral: 



T{gg) = u^{ia^vq^) m^AM + miBN + m^BH + m^iDS 



(A.24) 



where 



6(r,-l)2 



3ri-l 



3(r,--l)2 " 3(r,--l)2 ' 6(r,--l)2 
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2(ri-l) 



2(r,-l) 



H=-3l2 + l^^h - \-^h, 



s = 



2ri 



r/2- 



1 



3(r,-l)2" 3(r;-l)2" 3(r;-l)2 
For (A.20), note that: 

ka^ap^^ -kn{k + q)ii+gii^{2q + k)-k 
we simplify the structure: 

^i"+ii"+iii". 

where 



(A.25) 



(A.26) 



(A.27) 



ni"= (vtp,^vtp,),,(2..,^.^)l^±^^ 

Using the same method in diagram b we get the contribution of the above three terms: 
Tv = (io^yq^) [ \m];m,_iA^ ^ ^ 



(A.28) 



12(/:2-m2)2(/,2_^;„2^)(^2_^2^)3 

k"^ lis. 



1 



6(F-m?)3(F-^m2,)(^2_^2^)2j 4 ^ (^2 _ ^2)2(^2 _ ^^2^) (^2 _ ^2^)2 
+ ^^m, m,D _ _ ^;„2^) (^2 _ ^2^)3 / • 

Tir = (^W)-^-.^(-^)/(,2_^2)2(,24m2,)(^2_^2^)2- 

Till' = (jo-juv?^){m^A 



2(^2_^2)(^2_^^2^)(^2_^2^)3 



+ 



Thus we can get: 
T 



1 

2 (^2 _ ^2) (^2 _ ^^2^) (^2 _ ^2^)2 

1 ^2 

2 (^2 _ ^2) (^2 _ ^^2^) (^2 _ ^2^)2 

(k) 



+ miB 



rl 



]}• 



2(^2_^2)(^2_^^2^)(^2_^2^)3 

(A.29) 



iifi{i<yfivq^) m^AM' + miBN' + m^DW' 



Ufxi 



(A.30) 
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where 



(r,-l)(l-^)2 ^ 12(1-03 



-'3 H TTTT TTTT-, FTl M 



12(r,- 1)2(1-^3 



r/7 



/8 



6(r,-- 1)2(1-^2 I2(r,-- 1)2(1-0 6(1-^^ 



hi, 



4{n-\){\-^r ^'4(1-0^' 4(r,-l)(l-02" 2(r,-l)(l-0 



-rK4-2r, + r,-^-30 r^03-2j) 
12(r,-l)2(.^-l)2 i2(^_i)3^3 + 



/4- 



4(r,-l)(^-l)2 
/8. 



12(r,-- 1)2(^-1) 



12(r,-l)2(.^-l)3 
Now we cope with the last term (A.21). Simplify first the combination 



(A.31) 



(A.32) 



where 



iH4) {vlPL + viP,) Y^';_;fV ^ +II(m,). 
Km) - [VL^R + VRi-Ljmi {k + p)^-mf ' 

We can deal with the above two terms using the same method in diagram b.Thus we get the 
contribution of them to the loop integral: 

T[I{mi)] = (ia^vq^) j^{mfmnj ^ 

1 



(A.33) 



12 (^2 _ ^2)2(^2 _ ^^2^)2(^2 _ ^2^)2 
k"^ 1 



12 (^2 _ ^2)2(^2 _ ^^2^)(A:2 - ^2^)3 6 {k^-mf)\k^ - ^ml){k^ - m\f. 

. 1 k^ 
+mfB- 



+miB 



A{k^-m]f{k^-^ml){k^-mlf 
1 k^ 

2(^2_^2)2(^2_^^2^)(^2_^2^)2 
k^ 1 



2(^2_^2)^^2_^^2^)2(^2_^2^)2 2 (^ - m2) (^2 _ ^^2^) (^2 _ ^2^)3 
k' 



{k^-mUk^-^ml){k^-mlY 
1 k'' 



+mf m^D 



12 (k^-mjnk^-Wwnk^-'ni)^ 
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1 

+- 



1 



+ 



2 {k^-mfnk^ - ^ml){k^-mlY 2 {k^-mf){k^ - ^ml,f{k^-mlY 

1 k^ 

2 {k^ - m?) {k^ -^ml){k^- m^f {k^ - mf){k^ -^m^){k^- m^f. 



At last we obtain the following contribution to the loop integral: 



ilniiOnvq^) m^AM" + miBN" + m^DW" u^, 

L ± — I, J L J 

where 

6(r,-l)(l-^) ^^12(1-^)3 3 



12(1 -^)3(r,-- 1)2 



-/4 + 



12(r, - 1)2(1 - 12(1-^)2^ 6(1-^)2 



6(r,- 1)2(1-^)2 
hu 



h 



N" = 



2{\-^r' 4(r,--l)(l-<^)^^ 2(1-^)2^^ 

-n-l^+An^-r}t ^2 1 , 

4(r,-l)(l-^)2 4 2(1-^)2 5 + 2(n-l)(l-^)^' 



W" = 



2(^-1)2 
+ 



6(^-l)(r,-l)2 4(r,-l)(^-l) 
/3 



+ 
+ 



rm-2) ^(-r, + 2^) 
12(^-1)3^ 4(^-1)2 
r,--2r2 + 2r3<^-r3^2 ^,- + 2^ -4r,-^ +r2^ 
12(r,- 1)2(^-1)3 + 4(r,--l)(^-l)2 



/4 + 



2(^-iy 



r/5 



1 



2(n-l)(^-l) l2(^-l)(n-l)2 



1 



12"^^ (^-1> 



r/9. 



(A.35) 



(A.36) 



(A.37) 



Using (A.24),(A.30) and (A.36),we finally obtain the contribution from diagram (a): 



T (a) = eglu^{ia^vq^ 



i 1 



{4n 



2Z2~r"'" 



m 



w 



m„A 



M{a)+P{a)\n[4 + Q{a)\n[t, 



+miB \M'{a) + P'{a) ln[r/] + Q'{a) ln[Q 
+m^B M"{a)+P"{a)ln[ri] 
+m^D Moia) + PD{a)\n[ri] + QD{a)\n[^]^ |m^, 
where M{a),P{a),Q{a),. . . are listed in Appendix H. 



(A.38) 
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Appendix B. The calculation of the contribution from Z diagram 

We start with diagram (b): We begin to calculate the amplitude of diagram Z whose Feynman 
diagrams are shown in Fig 2. 

^'^'= (k-pr-mi ^^-'^ 

=yair-A y^^u ypU^ - 1-0 V; .2 

Jjt-4-ma "^^-ma \ {k-pY-t,m^2 
= lst + 2nd, 

4^ - iyLyL)ap -2sld^p,clp = {VlVR)ap -l^S^^ (B.2) 

where 

lp-4-ma If-ma 

2nd= -{l-m-^) . ,^ ^ y^T^^^-^) <v 1^2 - (B.3) 

-ma If-ma [k-py-^m^ 

Simplify the term Ist.Note that the terms Juiln can be dropped because they have no 
contribution to the loop integral.Thus we can get: 

-%^+2^y^^ + 8m«^^ 

[(k-qy-ml)][k^-miy ^ ' ^ 

Substituting (B.4) into (B.l) we can get the contribution of the term 1st to the loop integral: 



r(lst) = ^u^iia^yq") [m^A'{FF)+maB'{GG)+m^D'{WW)]u^. (B.5) 



where 



^1 ^ , l-4r« 



3{ra-iy 3{ra-l) 3{ra-iy 3{ra-iy 
iGG) =^T,--^T2, 



where = m^/m^, and r„ are integral functions listed in Appendix. 

Now we proceed to the term 2nd.To pick up the revelant terms, again we will use: 

■P){k-q)^^ {io^vq^)]/)0 - ^ {iG^ivq^)m*^lc^, 

k0{k-p){k-p) ^ ^{io^vq^W- (B.7) 

We can also use the replacement ^ — )■ when ^ is on the far right.Thus we can get the 
contribution of the term 2nd to the loop integral: 

2 

r(2nd) = ^u^iia^yq'') [m^A'iFF)' + maB'{GGy + m^iD'{WWyW(B.S) 
4cw 
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(FF)' = 



1 



ra 



5 ra 



2(ra-l)2 3ra-lJ 



T2 - 2^aT3 + 



+ 



1 ra 



+ 



ra 



3r«-l 6(r«-l)2 



4 ra ,7 
3 ra - 1 3 
ra 



T4 



3(ra-l)^ 3 3 



—mlra^Tu - ^m^ra^^Tn, 
3 6 



T2 + i^T,+ 



1 3rct rg 

2 2 2(5;- 1) 



7i 



(WW)' = 



"2^^" ^^^^ + 6(7^^' + 



3(ra-l)2^2^ 2 



2(ra-l) 
1 



a 

B' = y(c^ c^ Pr + c^ c^ Pr) 

i 

^ ~^{^aiJ.^aiJ.^L'\' (^jif^an^R)- 
a 

At last we get the contribution of the diagram {b)z 
2 

rmz] = JJ^I (4^)2 '^M (^'t^Mvg^) {m^A' [MM{b) + PP{b) lnr« + QQ{b) In 
+m«5' [MM' {b) + PP' (Z>) In r« + Qg' (Z?) In <^ ] 
+m^D' [MMd {b) + PPz) (Z?) In ra + (^) In ^ ] } , 

where MM{b) , PP{b) , QQ{b) , ... are hsted in Appendix I. 
Now replacing Z by G°,we get: 



(^)go = 



Acwmi 



- f H 



where 
H 



N = 



-fna{ca^iPL + cl^PR) 
1 



1 



Simplify N we get: 



(B.9) 



(B.IO) 



(B.U) 



(B.12) 
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Substituting (B.13) into (B.l 1) and using the same method in the diagram b,we can get the 
contribution of the to the amplitude: 

2 

T[bGo] = ^^j^u^{io^yq''){m^,A'[MM{c)+PP{c)\nra + QQ{c)\n^] 
+maB' [MM' {c) +PP'{c)lnra + QQ'{c)ln^] 

+m^,D' [MMd{c) +PPDic) Inra + QQd{c) ln<^] jw^u, (B.14) 
where MM{c),PP{c), QQ{c),. . . are listed in Appendix I. 
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Appendix C. The loop functions appearing in Fig 1 are defined and the results are listed 

' J {27t)Hk^-mj){k^-^m^){k^-ml,) {4nf m^l - ^){ri- l){ri- ^) ' 
d^k 1 / 1 r.lnr. 



' J {27t)\k2-mj){k^-ml,f (4;r)2m2,(l-n)V l-nJ' 

J _ r_d^ 1 / 1 / 

^~'{27t)\k^-mff{k^-^ml,)~ i^^f^^win-^)^ n-^J' 

r 1 / 1 A Inr,- \ 

""'J {2Ttf{k''-m^f{k''-ml)~ {^Ttfml{ri-\)\ n-l)' 

I = / 1 - ' 1 (i I "'^"^^ 

h-- 



h = 

h = 



i27ty (k^-mjy{k^-ml) {47ty m^in - 1) 

1 / 1 / nln^, 

{27t)\k^-mj){k^-^mlf ~ i^^f'^win-^)^ ^ n-^ 
d^k 1 __ 1 1 
{2Kf {k^-^mif ~ "(4^2 
d'^k 1 _ /II 
{2n)' {k^-mjf ~ "(4^2;^' 
d'^k 1 /II 



n2':)™2 ' 



' J {27t)\k^-mjf{k^-Ww? m^rn^{n-^)^y ^ {n-^) J' 

d'^k 1 / / r, + 5(^ 2r; + ^, <^ 



/lO 



\ 2t, n-t, nJ 



{2n)\k^-m})\k^-^mlf (4^f ^^.(r, - 2^ r,- ^ 

J^fc 1 / /I J__^ln|' 

(2;r)4 (/t2-m?)^()t2-^m2,) ~ (4;r) Wl^-/- ^2 ^ 2^^ 



' {2n)\k^-nij)\k^-Wwf (4;r)2m6,(r,-^)H 2r,- r,-^''^ 
1 / 1 / r,- 2r,ln^ 



''~J {27t)\k^-m^){k^-^m^f ~ {47lf2mUn-^fy ^ ' ^-^ )' 
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Appendix D. The loop functions appearing in Fig 2 are defined and the results are listed. 



Ti 



d^k 1 / (<^-l)r«lnra-(ra-l)^ln^ 



{Inf {k^ - ml){k^ - Wz) (k^ - ml) {47cf m|( 1 - ^ (r« - 1 ) (r« - ^ 



J. ^ r _a*k 1 1 / ^ ra\nra \ 

^ 1 {2Kf{k^-ml){k^-mlf {AKfrnl{\-ra)\ l-ra^' 

T = f d'k 1 1 ^^"TA 

^ J {2Tl)\k^-mlf{k'^-^ml) {ATlfrnl{ra-^)\ ra-^J' 

J ^ f \ ^ i 1 Inrg \ 

^ J {iTtf {k^-mlf{k'^-ml) {ATtfml{ra-\)\ ra-l)' 

T = I ^ - ' ^ fl I '"^"^^ 

d^k 1 ill 



76 = 
Ti = 
Ts = 



d'^k 1 /II 



{2Kf{k'^-mlf {Anf2raml' 
d^k 1 __ 1 1 
(2;r)" (yt2-m2)3 " ~~{A%f2'^z 



f d'k 1 / 1 ^ (ra + ^)ln^ 



7^10 = 



(2;r)4()t2-m2)2^^2_^^2)2 " (47r)2m4(ra-0^^ 

{27t)\k2-ml)\k^-Wzf (4;r)2m|(ra-<^f V 2^ r^-.^ ""raJ' 

" i (2;r)'^(;fc2_^2)3^^2_^^2) (4;rfm|(ra-<^f ^2 2ra ra-^r 

J {27t)\k2-mlfik^-^mlf (4;rfm|(ra-<^)n 2ra Va-^ ^ )' 

i (2;r)4()t2-m2)()t2-^m2f {Anf 2m\{ra- ^f^ ^ ra-^ )' 
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Appendix E. The loop functions appearing in diagram(b) are listed. 



23 



^^^^ ^ -5 + ri + y-2ri^ ^ n -rf + ^+ Ir^ ^n^^ 



24(^-ir(r,-l) 



P{b) = n 



-4 + 3<§ 



, ^(-2 + r,- + ^-2r,^+2^2) ^3 



Qib) 
M'{b) = 



1 



24(r,- 03(^-1)2 8(r,-0' 



24(^-1)2 12(<^-l)(r,-0 8(r,-^r 



8(^-l)(r,-l) 4(n-^). 



n 



8(^-lf 8(r,-^) 
1 



.8(r,-l) 4(r,-^) 



+ 



■ + 



'-«(l-2^) 



16(r,--^)3 



(r,-02 L8(r,-1)2 8(r,-^) 
r,^(-2r2-2r,^ + ^2) 
24(r,--04 



12(^-l)(r,-l)2 24(r,--l)(r,-^)(^-l) 



r,(2r,-3r,^+^) r,^(2^-l) 



24(r,- 1)3(^-1)2 24(^-l)2(r,--02 



8(^-1) 8(r,-0, 



+ 



-2^ + 1 _^2r2 + 2r,-^-^2n 



24(^-1)2 2A{ri-^Y 



Calculation of the anomalous magnetic moment ofjl 

Appendix F. The loop functions appearing in diagram(c) are listed. 



M(c) 



-r? - +2^^ + 3r,<^2 _ 3^3 _ ^. ^3 ^^4 2^. + rf + ^+ 



P{c) = 



24(r,- 1)2(^-1)2 



Q{c) = -P{c) + 
M'{c) 
P'{c) 



,24(^-l)2(r-l)' 24(r,-^)2(^-l)2 
-r2 + 2r,-^+r2^-^2_2,.^2_2^3 

24(r,-04 
r,(3-r,-2^) 

24(^-l)2(r,-l)3' 

1 r -n^-n + l^ ^ r,-<^ 



r,-^L8(r,-l)(^-l) 4(r,-0 



■ + 



!2'(c) = 
Md{c) 
Pd{c) -- 



^(-2<g+3r,) 
4(r,-^)(r,-l)(^-l) ■ 8(^-l)(r,-^f 

8(^-l)(r,-lf 8(r,-^)3' 



1 r 4r, -6<^+3r,<^ r,<^ 



r,-^L24(r,-l)(^-l) 4(r,-^)J 
1 



+ 



48(r,-^)3(^~l) ^ 



2^ + l)r, r,-2 F {-4rf + rf ^ - 2^^ + 4rt^ + n^^) 



^-IL 4(r,-l)(r,-0 8{n-iy 



1 



^ (-4r,^ ^ + 2n^^ + ^^ + 3rf-3rf^+rf^^) 



Qd{c) = 



(^-l)n24(r,-l)2 



2A{n-^Y 



(r,--On 8(^-1)2 8(r,-^; 

24(^-l)2(r,-^)4 
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Appendix G. The loop functions appearing in diagram(d) are listed. 



25 



M{d) 

P{d) 

M'{d) 

P'id) 

Moid) 



1 r;(4r?-5r,<^-5<^^ 



24 



= -Q{d) = 



rHi-ri + 2^: 



4 ' 



4(r,-<^) 
-Q'id) 



2 ' 



Poid) = -Qoid) 
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Appendix H. The loop functions appearing in diagram(a) are listed. 



M(fl) = 



_ 8 - 35r; + 66r? - 41r? + 2rf 3r,- - 2 



24{ri-iy 

r; 



P{a) 



+ 



+ 



12(r,- 1)2(^-1) 
Br? 



A{^-\){n-^){ri-\)^ A{n-\) 



^rj-3^+2n^-2^^ 1+5^ 



24{n - 1)'(^ - 1) 24(^ - l){n -^f 24(n- -^f 



Q{a) = 
M'{a) -- 
P'{a) = 
Q'{a) = 
M"{a) : 
P"(a) = 
Mo^a) 



m-\) 24(r,--0 



1 ^n^^-^^ + r,-^^?>{3n-\) 



1 



4(^-l)(r,-^) ^ 4(r,-l) 



+ 



(l-r^n 4(r,--^)2^2(l-r,-) 
1 ^2(-3r, + 2^+r,0 
4 (^-l)2(r,-^)2 ' 

ll-13r, 
'4(r,--l)2' 
r,(-5 + 6r,-) 
■ 2(r,-l)3 ' 

3r;-l ^2^ + 5r,-l 



8(r,-l)2^ 6(n-l)3 4(^-l)(r,-^) 

-l + 7r,-5r;^ + ^-2r?^) , ^ 



+ 



■ + 



r? 



24(n-l)3 



2 



"4(n-l)3 (r,-l)4^4(^-l)(r,-l)(r,-0 4(r,-^)2(^-l) 



^-1 



12(^-l)(n-^)2 + 24(^-l)(r,--^)3 4(r,--02(^-l)2' 



24(r,--l)4 12(^-^)2 2A{n-^y 
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Appendix I. The loop functions appearing in diagram Z are listed. 



MM{b) = 7-3r« + 6r2^-28r3 ^ 30r3 - 22r2 ^ - r^^^ _ ^3 



PP{b) -- 
QQ{b) 



12(ra-l)3 12(r«-^)3 
rli-5 + Srg) rg^ (lArj - ISvg^ + 10^^) 

rg^{2Arl-25rg^ + lO^^) 
6{rg-^r 



2(ra-l)2 2(ra-^)2 ' 



QQ\b) 



2(r«-l)3 2(r«-03 
r«<^(-2r« + ^) 



MM (h\ 5 + lrg-l2rl + 6rl , 18^ -22r2^ + llr^.^^ _ ^3 

"^"^^^'^ = 6U-^' ^ ^ 

r«(-2r3+8r2^-15r«^2^6^3) 



QQoib) 
MM{c) 

PP(c) = -QQ{c 
MM\c) 



6{rg-^r 
rg{-29rl + l9rg^+4^^) 
12(r«-^)3 

rg^i24rl-25rg^ + l0^^) 
6{rg-^r 

rai^rg-^) 



PP\c) = -QQ'ic) = 

rg{nrl-l9rg^+S^^] 



MM Die) 



12(ra-^)3 

rg{-2rl + SrU-l5rg^^ + 6^') 



PPoic) = -QQd{c) ^, , 
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